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Motivated by the physics of mobile triplets in frustrated quantum magnets, the properties of a 
two dimensional model of bosons with correlated-hopping are investigated. A mean-field analysis 
reveals the presence of a pairing phase without single particle Bose-Einstein condensation (BEC) 
at low densities for sufficiently strong correlated-hopping, and of an Ising quantum phase transition 
towards a BEC phase at larger density. The physical arguments supporting the mean-field results 
and their implications for bosonic and quantum spin systems are discussed. 
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The models of interacting bosons (with or without dis- 
order) have been a subject of active research. They are 
studied for a variety of reasons, coming from different 
experimental systems, such as Josephson junction ar- 
rays [1], 4 He in porous media [2], disordered films with su- 
perconducting and insulating phases [3] , or more recently 
in the context of atoms trapped on an optical lattice [4] . 
The interplay of interaction, disorder and kinetic energy 
leads to the ground states that can be a superfluid, a Bose 
glass, a Mott insulator or a supersolid [5-11]. In the con- 
text of spin models too, the Schwinger boson mean-field 
theories provide a useful description of magnetism in the 
bosonic language [12-15]. 

Over the last decade, bosons have also been used 
in the context of quantum magnetism to describe the 
magnetization process of gapped systems with a singlet 
ground state such as spin ladders, the triplets induced by 
the magnetic being treated as hard-core bosons. These 
bosons may condense, leading to the ordering of the 
transverse component of the spins, but they might as well 
undergo a superfluid-insulator transition, leading to mag- 
netization plateaux [16]. For pure SU(2) interactions, 
and without disorder, the common belief is that the only 
alternative, not realized so far in quantum magnets, is a 
supersolid, i.e. a coexistence of these phases. 

In this paper, we propose that there is another possibil- 
ity, namely a pairing phase without single particle Bose 
condensation. Our starting point is the observation that 
the effective bosonic model of a frustrated quantum mag- 
net such as SrCu2(B03)2 [17] contains, in addition to the 
usual kinetic and potential terms, a correlated-hopping 
term where a boson can hop only if there is another boson 
nearby, and that this term can be the dominant source of 
kinetic energy in geometries such as the orthogonal dimer 
model realized in SrCu2(B03)2 [18, 19]. While the possi- 
bility of bound state formation was already pointed out 
in that context, the consequences of the presence of such 
a term on the phase diagram at finite densities have not 
been worked out yet. 

For clarity, we concentrate in this paper on a mini- 
mal version of the model, but we have checked that the 




FIG. 1: The diagonal bonds (— t') denote a correlated-hopping 
process, and the nearest neighboring bonds denote the single 
particle hopping (— t). 



conclusions apply to the more realistic model derived for 
SrCu2(B03)2 [20]. This model is defined on a square 
lattice by the Hamiltonian 
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where 6J, b r are boson operators and n r = b\b r . t and t' 
are the measures of single particle and correlated-hopping 
respectively. A hard core constraint that excludes multi- 
ple occupancy should in principle be included. However, 
we will concentrate on the low density limit, where this 
constraint is expected to be irrelevant. So in the follow- 
ing we work with regular (soft core) bosons. 

Since the correlated-hopping term in Eq. (1) is quartic 
in the single particle boson operators, the simplest thing 
to do is a mean-field theory. Since the system gains en- 
ergy through correlated-hopping by having two particles 
nearby, a natural choice for a mean-field is the pairing 
amplitude. The particle density and the kinetic ampli- 
tudes are the other choices for the mean-fields. In the 
following, we formulate a mean-field theory in terms of 
these order parameters defined by: 

A = (blbl ±s ) (pairing amplitude) 
K = (blb r ±s), k = (bl± s b r ±s'} (kinetic amplitudes) 
n = (b].b r ) (particle density) 

where 8 =/= 8', and 6,8' — x,y. The particle density is 
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taken to be uniform, and the kinetic amplitudes real. In 
principle we can allow for an internal phase in the pairing 
amplitude (a non-zero phase between x and y direction 
bonds) like in the mean-field theory of t-J model in the 
context of the high-T c cuprates. Here, wc take the inter- 
nal phase to be zero (the extended s-wave pairing). The 
corresponding mean-field Hamiltonian has the following 
form: 
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where £q, £k and Ak are given by 

£ = 8t'L[A 2 + k 2 + n/t'] 

£ k = -2{t + At'n){cosk x + cosk y ) 

—8t'n cos k x cos k y — (p + 8i'fi/) 
Ak = 4t'A(cosfcc + cos k y ) 

The Hamiltonian Hmf can easily be diagonalized us- 
ing Bogoliubov transformation for bosons. The canonical 
free energy density for Hmf is given as: 
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where i? k = y/ ^ — 4 A 2 , is the quasi-particle dispersion, 
P = ^/ksT, and A = fi + 8t'n' is the effective chemical 
potential. Re-defining the chemical potential in this way 
makes k' a redundant order parameter in the mean-field 
theory. Note that fi and k' appear with right combination 
to give A as the new chemical potential in £k- Hence, for 
a given n, / is purely a function of A, k and A. 

The self-consistent equations for the order parameters 
can be written 
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where n c is the condensate density (the occupancy of 
the zero energy mode if any). Since the model is two- 
dimensional, n c — for T > 0. At T = 0, n c may or 
may not be zero, and solutions must be searched with 
two strategies: assume n c = and solve these equations 
for the unknowns (k, A, A), or assume there is a zero 
energy mode (which fixes A), and solve for the unknowns 
(k, A, n c ). In that case, the wave- vector corresponding 
to the zero-energy mode (k = here) must be excluded 
from the sum. If several solutions are found for a given 
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FIG. 2: The mean-field quantum phase diagram, n* is the 
critical density at which, for a given a = t' /(t + t r ), the single 
particle condensate density n c becomes zero. This marks the 
onset of the pairing phase with gapped quasi-particle excita- 
tions. 



density, the one with lowest energy should be chosen. In 
practice, we only found one solution for a given density. 
In general, these equations are solved by simple iteration. 
Note however that, when n c = 0, Eqs. (4) and (5) can be 
still solved by iteration, but Eq. (6) needs to be solved 
for A by some numerical method at each step of the 
iteration. 

Let us first discuss the T = results. The most re- 
markable feature is that it turned out to be impossible to 
find a solution with a non-zero condensate at low enough 
density unless t' is very small. In other words, as soon 
as the correlated-hopping is not too small, there is no 
single-particle BEC at low density. The critical density 
n* below which this is the case is plotted as a function 
of a = t'/(t + t') in Fig. 2. This figure calls for some 
comments. First of all, t'/t need not to be large for the 
effect to be observable, which ensures the relevance of 
the present discussion for the quantum magnets such as 
SrCu2(B03)2- Besides, the critical value n* is quite small 
even when t' dominates, and our approximation to treat 
triplets as soft-core bosons is expected to be good in the 
whole range of Fig. 2. Finally, densities of a few percent 
are definitely accessible in the context of quantum mag- 
nets, the density being equivalent to the magnetization 
relative to the saturation value. 

Next we turn to the nature of this non-BEC phase. 
Clearly this cannot be a commensurate insulating phase 
of the type observed before since it occurs for a range of 
densities. In fact, its nature is best revealed by looking at 
the order parameters. While n c = in the pairing phase 
and Eg a p = in the single particle BEC phase, the mean- 
field solution for k and A is non-zero on both sides. In 
Fig. 3 the behavior of n c , k and A is shown as a function 
of n for a = 0.99. It is not surprising that A is non-zero 
in the single particle BEC phase. In fact, the single par- 
ticle BEC state means (6 T ) ^ 0, which further implies 




FIG. 3: The variation of the quasi-particle gap and the mean- 
field order parameters as a function of n for a fixed a. While 
n c 7^ only in the single particle BEC phase (n > 0.072), the 
pairing amplitude A / in both phases. 



FIG. 4: The behavior of different order parameters as a func- 
tion of a for n = 0.05. For a > 0.78, the condensate density 
(n c ) vanishes, and the system goes into the pairing ground 
state. 



that (fctfot) ~ (fot)2 ^ o. Hence A will always be non-zero 
in the single particle BEC phase. The correct measure 
of the existence of the pairing (independent of the con- 
tribution from the single particle BEC) is A — n c . We 
know from the calculation [see Fig. 4] that for the non- 
interacting Bose gas (a = 0), A = n c , as it should be. 
However, for any finite a we find A— n c > 0. Thus, for ar- 
bitrarily small values of the correlated-hopping, the sys- 
tem develops a tendency towards pair formation. How- 
ever, it does not suppress the single particle BEC in fa- 
vor of a purely pairing phase until sufficiently strong a 
is reached for sufficiently small n. 

The results of our mean-field calculation arc similar 
to those obtained on a different problem in the context 
of the atomic gases [21]. These are studies regarding 
the transition from a purely molecular condensate (MC) 
to an atomic condensate with a non-zero fraction of the 
molecular condensate present (AC+MC) accross the Fes- 
hbach resonance. Our pairing phase is like the MC, 
and the single particle BEC phase is analogous to the 
AC+MC. 

The nature of the transition based on symmetry con- 
siderations is also similar. The mean-field Hamiltonian, 
Hmf, explicitly breaks the U(l) gauge symmetry, how- 
ever it is still invariant under global Z 2 symmetry, that 
is under b r — > — b r . In other words gauge symmetry, 
b\ — > b\e % ^ leaves Hmf invariant for <f> — tt. This resid- 
ual Ising like symmetry will also be broken if there is 
single particle BEC (because (b) ^ 0). What we have 
in Fig. 2 is such an Ising symmetry breaking quantum 
phase transition, where n c is the relevant order parame- 
ter. The pairing phase respects this Z 2 symmetry while 
the single particle BEC phase breaks it spontaneously. 

The temperature dependence of various quantities in 
the mean- field theory is shown in Fig. 5. The tempera- 
ture at which A becomes zero is called T c . This quantifies 
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FIG. 5: The temperature dependence of A, k and n for a 
given fi and a. Inset: dependence of T c on the density. 



the mean-field phase transition from a normal Bose gas 
at high temperatures to a pairing phase below T c . The 
inset of Fig. 5 shows T c as a function of n. Remarkably, 
there is no detectable anomaly upon going through the 
critical density n*. Since we are in 2D, we do not expect 
to have a true BEC of pairs, but rather a Kosterlitz- 
Thouless (KT) transition. These results suggest that the 
system should undergo one KT transition whatever the 
density, followed by an Ising transition if n > n* . 

To check the validity of the mean-field approximation, 
hence of our conclusions, it would be very useful to have 
unbiased numerical results on the model of Eq. (1). How- 
ever, we have good reasons to believe that the predic- 
tions of the present mean-field theory are physically rel- 
evant. Mathematically, the structure of the mean- field 
equations and the results of the calculations are similar 
to the Schwinger boson mean-field theory of the quantum 
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spin system [13, 15]. In that context, the single particle 
BEC phase implies an ordered phase in the spin variables, 
while the pairing phase denotes a disordered phase. Now 
the physical relevance of these disordered phases is well 
established in the context of quantum magnets [22, 23], 
and we expect the same to be true here. 

The model studied in the present paper has similari- 
ties with the ring exchange model of bosonic Cooper pairs 
introduced some time ago by Paramekanti and collabo- 
rators [24]. In their model, the pairs of bosons hop on 
the opposite corners of a plaquette. The low temperature 
physics is significantly different however. In their model, 
the number of bosons is conserved on each row and col- 
umn of the square lattice, leading to the Luttinger liquid 
like physics and critical correlations in the ground state. 
In our model, the correlated-hopping does not sustain 
any such conservation law, and the ground state is ex- 
pected to develop a true long range order. 

Finally, let us briefly discuss the physical implications 
of these results for the magnetization process of gapped 
quantum magnets. The thermodynamics was already 
discussed in the boson language: we expect to observe 
a KT transition for any magnetization, followed by an 
Ising transition if the magnetization is larger than a crit- 
ical value. This will remain essentially true for 3D sys- 
tems, the KT transition being replaced by a true phase 
transition toward an ordered phase. However, we also ex- 
pect very significant differences between the zero temper- 
ature phases. Single particle BEC means magnetic long- 
range order, and the system is expected to have gapless 
transverse spin waves. In other words, the gap detected 
in spectroscopies such as inelastic neutron scattering or 
NMR will vanish. However, at low magnetization, we 
only have pair BEC. The order implied by this pair BEC 
will be of nematic type since the transverse components 
of the spins within a pair can be flipped without chang- 
ing the correlations. But more importantly, there is a gap 
to single particle excitations, i.e. to single spin flips. Al- 
though the system is gapless in this phase, we thus expect 
to observe a gap in neutron scattering or NMR, the gap- 
less excitations appearing only in the channel AS = 2. 

In summary, we have shown that the correlated- 
hopping can change drastically the properties of bosons, 
leading at low densities to a pairing phase without single- 
particle BEC, and with gapped quasi-particle excita- 
tions. In the context of quantum frustrated magnets, this 
leads to the prediction of an Ising phase transition (for 
low magnetization) as a function of the magnetic field, 
for systems where frustration reduces direct hopping of 
triplets, thus making correlated-hopping the dominant 
process of kinetic energy. 

Beyond frustrated magnets, these results will have im- 
plications on all systems where correlated-hopping may 
be the dominant source of kinetic energy. One such class 
of systems are the atomic gases, where different external 
parameters control the hopping and the Coulomb terms. 



Whether instabilities of the kind described here can be 
induced in these systems by reducing the single-particle 
kinetic energy is left for the future investigation. 
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